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Abstract 

The QCD corrections to the moments of the invariant mass distribution 
in the semileptonic r decays are considered. The effect of the renormalization 
scheme dependence on the fitted values of a s (m%) and the condensates is dis- 
cussed, using a simplified approach where the nonperturbative contributions 
are approximated by the dimension six condensates. The fits in the vector 
and the axial-vector channel are investigated in the next-to-leading and the 
next-to-next-to- leading order. The next-to-next-to- leading order results are 
found to be relatively stable with respect to change of the renormalization 
scheme. A change from the MS scheme to the minimal sensitivity scheme 
results in the reduction of the extracted value of a s (m^.) by 0.01. 
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1 Introduction 



Recently there has been considerable progress in the determination of the vector 
and axial-vector hadronic spectral functions from the semileptonic decays of the 
t lepton Using the QCD predictions for the moments of these spectral 

functions one may obtain constraints on a s and the parameters characterizing the 
nonperturbative QCD dynamics (|0]-[0, ^B|-|]T3|). The accuracy of the obtained 



value of a s appears to be quite high. In order to have a proper understanding of 
the phenomenological relevance of this determination of a s it is important to make 
a careful estimate the theoretical uncertainties in the QCD predictions |I6|-|]26|. 
In this note we investigate the uncertainties in the evaluation of the perturbative 
part of the QCD prediction for the spectral moments, for which the next-to-next- 
to-leading order (NNLO) approximation is available. We study in some detail the 
sensitivity of the perturbative QCD predictions to the choice of the renormalization 
scheme (RS) and the effect of the RS dependence on the fitted values of the QCD 
parameters. 



2 Theoretical framework 

Let us begin by summarizing the theoretical framework adopted in the analysis of 

V/A 



the r decay data |j7f]-[]12]. In it was suggested to use the R^v/a moments, defined 



by the relation 



' T > v/A T P Jo I ml [ml ds 




where dT^ d A /ds denotes the invariant mass distribution for the Cabbibo allowed 
semileptonic r decays in the vector (V) or axial-vector (A) channel and T e denotes 
the electronic width of the r lepton. The QCD predictions for R^y/A nave the form: 

3 

R t,v/a = 2 \ Vud \ Sew R (1 + $pt + S npttV / A ), (2) 

where V u d is the CKM matrix element (\V u d\ = 0.9752) and Sew = 1-0194 represents 
the correction from electroweak interactions ]^7[ P^] . The Rq factor denotes the 
parton model prediction — for later use we shall need Rq° = 1 and R^ 2 = 13/210. 
The 5pl term denotes the perturbative contribution, evaluated for 3 massless quarks. 
(The u and d quark mass effects are negligible, the s quark mass effects are very 
small for non-strange decays and the c quark is considered to be decoupled [^9], |3~0~||.) 



The 5$. term is universal for the V and A channels. The 8^ pt V j A term in (|2|) denotes 
the contribution from the nonperturbative QCD effects, which are estimated using 
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the SVZ approach 



$npt,V/A - Y ^{D)V/A— Y ZTd X] C D,j < >, (3) 

D=4,6... D=4,6... t j 

where < Opf > are the vacuum expectation values of the gauge invariant operators 
of dimension D and • are coefficients specific for the considered spectral moment 
and the type of the operator. The • coefficients are in principle power series in 
the strong coupling constant. 

The object of greatest interest is of course the total decay rate for the Cabbibo 
allowed semileptonic r decays in the vector/axial- vector channels, R®° v / A . The per- 
turbative correction to this moment is sizeable and it is highly sensitive to the value 
of the strong coupling, due to the low characteristic energy scale of m T = 1.777 GeV. 
The higher R* 1 moments are introduced to take advantage of the the full information 
contained in the hadronic spectral functions. By using the predictions for several 
jjki moments one may obtain a simultaneous fit of a s (m^.) and of some of the con- 
densates < Op f >. In this way the whole analysis becomes self-consistent and 
in addition one obtains a check on both the perturbative and the nonperturbative 
QCD contributions. 

The perturbative QCD corrections 5^ are evaluated using a contour integral 
expression which relates them to the QCD correction 5^ to the so called 



Adler function [|34|, i.e. the logarithmic derivative of the transverse part of the 

r (i)i 



vector/axial- vector current correlator H^ 1 '• 



-12 7 rV^ni 1 j^(a) = 3[l + <(-a)]. (4) 



(In the approximation of massless quarks the perturbative contributions to the Adler 
functions for the vector and axial- vector current correlators are identical.) We have: 

« = - I -f kl ( 4l *£{-*), (5) 



pt tc Jc <j \ m% j 

where f kl (a/m^.) is a weight function specific to the considered moment and C is a 
contour running clockwise from a = m 2 ^ — ie to a = + ie away from the region 
of small |er|. In the following we shall need the weight functions / 00 and f 12 : 

f\x)= l --x + x*- l -x\ (6) 

*i2, \ 1 70 3 105 4 126 5 175 6 60 7 

f u (x) = x j H x 4 H x 5 x 6 H x 7 . (7) 

J K 1 2 13 26 13 13 13 y ' 
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The NNLO renormalization group improved perturbative expansion for 8^ may 
be written in the form: 

a£(-<r) = a(-a)[l + na(-a) + r 2 a 2 (-a)}, (8) 

where a = a s /n = g 2 /(47r 2 ) denotes the running coupling constant that satisfies the 
NNLO renormalization group (RG) equation: 

a— — —- a { 1 + cia + c 2 a ). (9) 
acr 2 

In the modified minimal subtraction MS scheme (i.e. using MS subtraction pro- 
cedure and choosing fi 2 = —a) we have []3~5|-||38|| for n/ = 3 rf IS = 1.63982 and 



r^ s = 6.37101. The renormalization group coefficients for rij = 3 are b = 4.5, 
ci = 16/9 and S* 3 = 3863/864 « 4.471. 

The QCD predictions for the S^j. are usually calculated in the modified minimal 
subtraction (MS) renormalization scheme |39|| . However, in the NNLO approxima- 
tion with massless quarks there is a two-parameter freedom in choosing the RS. 
This is a consequence of the fact that in each order of perturbation expansion the 
finite parts of the renormalization constant for the coupling constant may be cho- 
sen arbitrarily. Different choices of the finite parts of the renormalization constant 
result in different definitions of the coupling constant, which are related by a finite 
renormalization. This results in a change of values of the coefficients r 1; r 2 and c 2 . 
(We restrict our discussion to the class of mass and gauge independent schemes, 
for which the coefficients b and c\ are universal.) The formulas describing how the 
redefinition of the coupling affects the coefficients r« and c 2 are collected for ex- 



ample in [40]. The dimensional QCD parameter A also depends on the choice of 



the RS 41 . In the NNLO there exists however a RS invariant combination of the 



expansion coefficients | 4"2"|| 



p 2 = c 2 + r 2 - c\r x - r x . (10) 

For the 5n we have p 2 = 5.23783. 

The change in the expansion coefficients and the change in the coupling constant 
compensate each other, but of course in the finite order of perturbation expansion 
such compensation may only be approximate, which results in the numerical RS 
dependence of the perturbative predictions. This RS dependence is formally of 
higher order in the coupling constant, but in the case of the r decay it may be 
significant numerically, since the coupling is not very small at the energy scale of 
m T . It is therefore very important to verify to what extent the RS dependence 
affects the predictions and the fits to the experimental data. 

The authors of [[I], 0, £| used the Rf moments with k = 1 and / = 0,1,2,3, 
and fitted the D = 4, 6, 8 condensates. Since the aim of this note is primarily to 
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study the theoretical uncertainties in the whole procedure we shall adopt a simplified 
approach, which still has considerable phenomenological relevance. This approach 
is based on the fact that the dominant nonperturbative contribution to R^°v/a in 
the SVZ expansion comes from the D = 6 term |7], |J, because the D = 4 term 
is suppressed by additional power of a s . We shall therefore neglect the D = 4 
contribution to the total decay rate and - - for consistency - - the higher order 
correction to the D = 6 coefficient. (Such an approximation was in fact made 
already in In the following we shall also neglect contributions from the D > 8 
condensates. In order to be able to fit the D = 6 condensate together with ot s 
we shall use the QCD prediction for the -R^V/a moment, which similarly to -R°°y/A 
has a suppressed contribution from the D = 4 condensate. Neglecting in R^ v /a 
the contributions from D > 8 condensates we obtain a simple set of self-consistent 
formulas. In our approximation: 



A 00 _ A 00 _ 2471 " 2 \ " r ^ n V/A {u , 

°npt,V/A - °(6)V/A- — e~ 2^ G 6 ,j < ^e,j >> 1 11 J 

m T j 



*12 _ Jfl2 16807T^ ^ v/A no x 

°npt,V/A - °(6)V/A - 13 m6 °6J < <^6,j >> U^J 



where X)j Cej < ^6j A > i s the leading D = 6 contribution to the transverse part of 
the hadronic vacuum polarization function: 



-v) 3 ^U*) = E^<0^>. (13) 



This contribution is dominated by the four-quark condensates flO |. In the phe- 
nomenological analysis it is usually expressed in a simplified form, motivated by the 
chiral symmetry and the vacuum saturation approximation [|I~C|j : 



C 6 j < V Jf >= h v/A ^-a s p < qq > 2 , (14) 

3 

where < qq > is the quark condensate, hy = —7, — 11 and p is an effective 
parameter, characterizing the deviation from the strict vacuum saturation. 

In our study of the RS dependence effects we parametrize the freedom of choice 
of the RS by the parameters T\ and C2, following the conventions of our previous 
work E4l p5]. To obtain the perturbative QCD corrections we evaluate the 



J pt 

contour integral numerically, using under the integral a numerical solution of the 
RG equation @ in the complex energy plane. In this way we take full advan- 
tage of the RG-invariance properties of the perturbative prediction and we resum 
to all orders some of the large terms which would otherwise appear in the pertur- 



bation expansion [46, 47]. We assume that the integration contour C is a circle 
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(j = —m 2 exp(—i6), 9 G [— it, 7r]. To determine the numerical value of the running 
coupling constant on the contour C we solve the transcendental equation, which 
results from integration of the RG equation (§J) with a suitable boundary condition 
and analytic continuation to the complex energy plane: 

bin (^) -i b l = r f^-r 1 + c 1 In (±\ + F<»>(a), (15) 

V MS J ^ X ' 

where in NLO 

F^(a) = - + Cl \n(-^-), (16) 
a Vl + Cia/ 

and in NNLO for 4c 2 - c\ > 



1 c 
F (2) (a) = - + ci ln(cia) - — ln(l + c x a + c 2 a 2 ) + 

"^-^ arctanf"^-^' 72 ). (17) 



(4c2-c?)V2 V 2 + Cl a 

The presence of A-^ in the expression valid for arbitrary scheme follows from our 
taking into account explicitly the relation between A parameters in different schemes 



HH and usine; A-W^ reference parameter. 

After evaluating the predictions for 5^ and b p l we perform the fits to the experi- 
mental data for R°^ v / A and D^ V , A = R^v/a/ Rt°v/a- We use the experimental values 
reported recently by ALEPH g: R™ v = 1.782 ± 0.018, D™ v = 0.0532 ± 0.0007, 
roo a = 1.7H ± 0.019, D\ 2 A = 0.0639 ± 0.0005. In the fits we' assume for simplicity 
that the experimental errors for R®° and D^ 2 are not correlated. 

We express the results of our fits in terms of a s (m 2 ) and S^y (We actually fit the 

value of A|^, which we then convert to a s (m 2 ) using the RG equation in the MS 
scheme.) For comparison with other determinations of the strong coupling constant 
we extrapolate a s from m 2 to m 2 z . Our procedure for extrapolation relies on the 
matching formula relating a s (fj J 2 ,nf + 1) to a s (/i 2 ,ny): 

/ 2 t \ / i \ LaK^^f) 1 / - 57 _ 11\ al (/i 2 , n f ) 

« s (/i 2 , n/ + 1) = a s (fi 2 , n/) + - ^ /7 +-[L 2 + -L- —J /; , (18) 

where L = ln(/i/m g ) and m 9 is the running quark mass m q (fi 2 ) of the heavy quark 
evaluated at the scale fi = m q . (The NNLO matching formula of that form was 



originally proposed in |48|, |||, ^(J, see also discussion in [pi]] . However, in |30| it was 
found that a numerical coefficient in the NNLO term is actually different, which was 
subsequently confirmed in | 52|] . We use the coefficient of @.) In order to evolve a s 



from scale fii to the scale /z 2 we solve the equation: 

F^(a s (fil)/7r) - F^(a s (fil)/n) = Mn^/^). (19) 
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Table 1: Table of values of a s (m 2 ) and a s {m 2 z ) in NNLO related by the matching 
procedure described in the text. 

To obtain a s (m 2 z ) from the given value of a s (m 2 ) we first evolve a s from the scale 
m T to the scale of 2m c using the rif = 3 RG equation, then we use the matching 
formula (|l"8"D to obtain a s ((2m c ) 2 ,n; = 4), evolve this to the scale of 2m b using 
the 7if = 4 RG equation, use the matching formula to obtain a s ((2m b ) 2 , rif = 5), 
and finally evolve this to the scale of mz using the n/ = 5 RG equation. We use 
rh c — 1.3 GeV and = 4.3 GeV, which are the central values recommended by the 
Review of Particle Properties [53]. In Table [l] we give for reference some values of 



a s {m 2 ) and a s {m 2 z ) as a function of m(m T /A^y in the NNLO approximation. 



3 Fits in the vector channel 

Let us begin with the calculation in the MS scheme, to see how our approximate 
treatment compares with a more complete analysis reported in [f|]. In Table ^| we 
give the values for SB in the MS scheme, in NLO and NNLO, as a function of 

ln(m T /A^-). For completeness, we also include precise values for 5°°, which may 



be compared with those given previously in |10| . Fitting the experimental results 



for R° T ° V and D X2 y we obtain in NNLO = 441 ± 32 MeV, which corresponds to 
a s {m 2 ) = 0.367±0.018 and a s (m 2 z ) = 0.1238±0.0018. This is very close to the value 
a s (m 2 ) = 0.360 ± 0.022 obtained by ALEPH in a fit involving more Rf moments 
and D = 4, 6,8 nonperturbative contributions [|J. We also obtain S^ v = 0.0147 ± 
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Table 2: Table of values of the NLO and NNLO predictions for 5®® and 5^ in the 
MS scheme. 



0.0025, which is in reasonable agreement with the value of the nonperturbative 
contribution obtained by ALEPH Q . This confirms our expectation that the D = 6 
approximation adopted in this work provides a good approximation to the more 
complete analysis involving larger set of parameters. 

Let us now consider the same fit, but in a different renormalization scheme. 
As is well known, the theoretical and phenomenological motivation for the widely 
used MS scheme is not very strong, and there has been extensive discussion on the 
problem of the optimal choice of the renormalization scheme |4"2"| , |43| , [P], [45] , [[54|]- 



jc>8|j . One of the interesting approaches is based on the so called principle of minimal 
sensitivity (PMS) [fl"2"| |. The philosophy behind this approach is very simple — since 
the theoretical predictions of any theory should be in principle independent of the 
RS, then in the finite order of perturbation expansion one should look for the RS, 
which mimics this as close as possible. 

In the case of the conventional perturbative QCD expansion the RS parameters 
of the PMS scheme are determined by a system of transcendental and algebraic 



equations |42j| . Unfortunately, in the case of perturbative predictions obtained via 
numerical evaluation of the contour integral these equations do not apply, so the 
optimized parameters have to be determined by direct computation of 5 kl for differnt 
values of r% and c 2 . 

The dependence of <5°° on the scheme parameters r\ and c 2 was discussed in 



detail in |£J and the RS dependence of tffi* was investigated in |23], |25|. In both 
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Table 3: Table of values of the NLO and NNLO predictions for obtained in a 
scheme preferred by the principle of minimal sensitivity (r x = —0.64 in NLO and 
ri = 0, c 2 = 1.5p 2 = 7.857 in NNLO). 

(3) 

cases it was found that for moderate values of A12^ the NNLO predictions have a 
saddle point type of behavior as a function of r± and c 2 and that the position of the 
saddle point is well approximated by r\ = and c 2 = 1.5p 2 = 7.857. (Incidentally 



these scheme parameters correspond to the approximate solution |60[ of the algebraic 



PMS equations for 5n evaluated for spacelike momenta.) For very large values of 

(3) 

A.)-f- the RS-dependence pattern is more complicated than a simple saddle point, 

Ivl O 

but even then the scheme parameters distinguished above belong to the region of 
extremely small RS dependence. We shall therefore accept these parameters as the 
PMS parameters in NNLO. The values of the NLO and NNLO PMS predictions for 
5^ are given in Table || (The values for 5®® have been already given in [p3|| . The 
contour plots provided there in principle allow one to obtain predictions for 5®® in 
arbitrary scheme with reasonably large expansion coefficients.) 

Using the PMS predictions we obtain from the NNLO fit in the vector channel 
<5°° )vr = 0.0156 ± 0.0023 and A§L = 421 ± 30MeV, which corresponds to a s (m 2 T ) = 
0.356 ± 0.017 and a s {m 2 z ) = 0.1226 ± 0.0018. We see that in NNLO the change 
from the MS scheme to the PMS scheme results in the reduction of the fitted value 
of a s by an amount significant compared for example to the presently available 
experimental precision. 

In order to make our calculations more generally useful we show in Fig. |I] the 
results of the NNLO fit of a s (m 2 ) and 6^ v , obtained using the PMS predictions, as 
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a function of the experimental values of -R°°y and D^ v . For comparison we indicate 
the results of the NNLO fit in the MS scheme (dashed lines). 

. 4 1 




0.012 0.014 0.016 0.018 0.02 0.022 0.024 



Figure 1: Plot of the fitted values of a s (m^) and 5^ in the vector channel as a 
function of R®° v and -D*y, obtained using the NNLO PMS predictions. The dashed 
lines indicate the change in the plot when the MS NNLO predictions are used instead. 

Besides PMS another frequently discussed approach to the optimization is the 
effective charge (EC) method f£5|| , which amounts to the absorption of all the 
higher order radiative corrections to the physical quantity — in this case <5n — into 
the definition of the coupling constant. In the EC scheme we have ri = and 
C2 = p2- The NLO and NNLO predictions for 5^ and 5^ in the EC scheme are 
given in Table [|. We see that in NNLO the difference between EC and PMS is very 
small. This is reflected by the results of the fit: in NNLO we obtain a s (m^.) = 0.356 
and 6$ )v = 0.0158. 

To have a broader picture of the renormalization scheme dependence we perform 
the fit in a more general class of schemes. It is clear that regardless of our choice 
of the concrete optimal scheme there is a continuum of schemes, which seem to be 
equally reasonable. Predictions in such schemes should also be somehow taken into 
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Table 4: Table of values of the NLO and NNLO predictions for 5®® and in the 
Effective Charge scheme. 

account in the phenomenological analysis. A natural way to do this is to supple- 
ment the prediction in a preferred scheme with an estimate of the variation of the 
predictions over a whole set of a priori acceptable schemes. A concrete realization 
of this idea was presented in |59| , based on the existence of the RS invariant p 2 , 



which provides a natural RS independent characterization of the magnitude of the 
NNLO corrections for the considered physical quantity. In |5{| it was proposed to 
calculate variation of the predictions over the set of schemes for which the expansion 
coefficients satisfy the condition: 

o-2(n,r 2 , c 2 ) <l\p2\, (20) 

where 

o-2(ri,r 2 ,c 2 ) = |c 2 | + |r 2 | +c 1 |r 1 | + r\. (21) 

A motivation for the condition ( p0|) is that it eliminates schemes in which the ex- 
pansions (H) and ([D involve unnaturally large expansion coefficients that introduce 
large cancellations in the expression for the RS invariant p 2 . The constant I in 
the condition (|20|) controls the degree of cancellation that we want to allow in the 
expression for p 2 . In our case we have for the PMS parameters cr 2 (PMS) ~ 2|p 2 |, 
so in order to take into account the schemes, which have the same — or smaller 
- degree of cancellation as the PMS scheme we take 1 = 2. One may expect, that 
the estimate of the RS dependence obtained according to this prescription would 
be useful for a quantitative comparison of reliability of perturbative predictions for 
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different physical quantities, evaluated at different energies. It is also clear that any 
large variation of the predictions over a set of schemes satisfying the constraint (|20| ) 
with I = 2 would be an unambiguous sign of a limited applicability of the NNLO 
expression. 

In Fig. H] we show how the value of a s {m^) resulting from the fit depends on the 
parameters r\ and C2 specifying the renormalization scheme in NNLO. In the region 
of scheme parameters satisfying condition (^) with I = 2 the minimum value of 
a s {ml) = 0.347 (A§L = 403 MeV, a s (ml) = 0.1217) is attained for n = -1.62 and 

1 VI o 

C2 = 0, and the maximum value of a s {m 2 .) = 0.367 (A^- = 440 MeV, a s (m 2 z ) = 
0.1238) is attained for r± = 0.96 and C2 = 0. It should be noted that the MS scheme 
parameters lie outside the / = 2 region, but the fitted value of a s (m^) in this scheme 
coincides with the maximal value quoted above. It is also important to note that 
one of the other potentially interesting schemes, for which the NNLO expressions 
recently became available |)2] — the so called V scheme |61j] - corresponds to the 



RS parameters r\ = —0.109 and C2 = 26.200, which lie very far outside the I = 2 
region. It seems that the large value of C2 in the V scheme restricts its usefulness 
for the low energy perturbative predictions, because of the influence of the Landau 
pole in the RG equation. 

It is of some interest to perform the same fits using instead the NLO predictions. 
Using the NLO predictions in the MS scheme, we obtain 5^ v = 0.0148 and A|^ = 
527 MeV, which translates via NLO RG equation into a s (m 2 ) = 0.388, which in 
turn corresponds via NLO extrapolation to a s (m 2 z ) = 0.1261. We see that the 
value of is practically identical to that obtained in the NNLO fit. However, 
the value of the strong coupling constant is surprisingly high and it is significantly 
different from the NNLO value. Using the NLO predictions in the PMS scheme we 
obtain = 0.0150 and a|^- = 465 MeV, which corresponds to a s {m 2 ) = 0.358 

and a s (m z ) = 0.1230. We see that although the difference in the NLO and NNLO 

(3) 

values of A^- obtained in the PMS fits is still considerable, it is nevertheless much 

smaller than in the case of the MS scheme. (If we look at a s (m T ) instead of A-^ 
this difference is even smaller, although this may be partially a result of a fortuitous 
compensation between different factors influencing the evaluation of a s in various 
orders.) For completeness we give the NLO results in the EC scheme: S9^, v = 0.0149 

and AgL = 472 MeV (a s (m^) = 0.361). These numbers provide a nice illustration 
of the fact, that the difference between the NLO and NNLO predictions is strongly 
scheme dependent. Therefore, if such a difference is to be used in any way to estimate 
the precision of the QCD prediction, some way of making a preferred choice of the 
renormalization scheme must be employed. 

So far our discussion was concentrated primarily on the value of a s {m 2 .) coming 
from the fit. In Fig. |3| we show the RS dependence of the fitted value of S9^ v . We find 
that in the set of schemes satisfying the constraint ( p0|) with / = 2 the fitted value 
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r i 

Figure 2: The contour plot of the fitted value of a s (m 2 T ) in the vector channel as a 
function of the RS parameters r\ and c-z- For technical reasons we use C2 — c\T\ as 
an independent variable on the vertical axis. The region of the scheme parameters 
satisfying the condition (p0|) with I = 2 has been also indicated. 



of 5®Qj V changes in the range 0.0145-0.0182. Using the Eq. (|14"D we may translate 
this into the range (2.22 — 2.78) x 10~ 4 GeV 6 for the commonly used parameter 
a sP < QQ > 2 - This seems to be in reasonable agreement with the values obtained 
previously by other authors, for example 1.8 x 10~ 4 GeV 6 obtained in the original 
work of SVZ PJ or (3.8 ± 2.0) x 10" 4 GeV 6 obtained in a more recent analysis |63| 



4 Fits in the axial-vector channel 

Similarly as in the case of the vector channel we start with the fit in the MS scheme. 
We obtain S$ )A = -0.0168 ±0.0021 and A§L = 398±37MeV, which corresponds to 
a s (m 2 T ) = 0.344 ±0.019 and a s (m|) = 0.1213 ± 0.0021. The result of this fit should 
be compared with ot s (m 2 ) = 0.365 ± 0.025, obtained by ALEPH in a more complete 
fit [|]]. We find that the difference between our results and the ALEPH result is 
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Figure 3: The contour plot of the fitted value of 5^ in the vector channel as a 
function of the RS parameters T\ and ci- The region of the scheme parameters 
satisfying the condition (B^) with I = 2 has been also indicated. 



slightly bigger than in the case of vector channel. It should be noted however that 
the ALEPH fit in the axial-vector channel has surprisingly large D = 4 contribution, 
which would be difficult to justify theoretically. 

Performing the NNLO fit in the PMS scheme we obtain S^ A = —0.0165 ± 

0.0018 and = 380 ± 34 MeV, which corresponds to a s (ml) = 0.335 ± 0.018 and 
a s (m z ) = 0.1203 ± 0.0021. Similarly as in the case of the vector channel we find 
appreciable reduction in the extracted value of a s (m^) as compared to the values 
obtained in the MS scheme. 

To make our results useful in the case of future improvements of the experimental 
analysis we show in Fig. |] the results of the NNLO fit of a s (m^) and S^ A , obtained 
with the PMS predictions, as a function of the experimental values of R^ A and D^ 2 A . 
For comparison we indicate the results of the NNLO fit in the MS scheme (dashed 
lines) . 

In Fig. |5] we show how the fitted value of a s {m 2 ) depends on the RS parameters 
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Figure 4: Plot of the fitted values of a s (m^) and 5^ in the axial- vector channel 
as a function of R®°a an d D^ A , obtained using the NNLO PMS predictions. The 
dashed lines indicate the change in the plot when the MS NNLO predictions are 
used instead. 



T\ and C2- In the region of the scheme parameters satisfying the condition ( ^0|) with 
I = 2 we have 0.326 < a s (m 2 T ) < 0.343 (364 MeV < < 397 MeV, 0.1193 < 

a s {m 2 z ) < 0.1212). Similar figure may be obtained for the RS dependence of 69q\ a : 
it is found that 0.015 < -S$ )A < 0.017. 

Having obtained the results for 5?^ in the vector and axial-vector channels it is of 
some interest to verify the simple relation between them, implied by the generalized 
vacuum saturation approximation: 5^ A /5^ V = —11/7 ~ —1.57. Using the num- 
bers from the NNLO PMS fits in both channels we obtain S^ )A /6^ )V = -1.06±0.28. 

In order to have a full picture of the perturbative uncertainties in the axial- 
vector channel we also perform the NLO fits. Using the MS scheme, we obtain 
$(6)A = — 0.0166 and = 473 MeV, which translates via NLO RG equation into 
a s {m 2 ) = 0.362, which in turn corresponds via NLO extrapolation to a s {m 2 z ) = 
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Figure 5: The contour plot of the fitted value of a s (m%.) in the axial- vector channel 
as a function of the RS parameters r% and C2. The region of the scheme parameters 
satisfying the condition (^) with / = 2 has been also indicated. 

0.1235. In the PMS scheme we obtain 5g? )A = -0.0167 and = 419 MeV, which 
implies a s (m^) = 0.336 and a s (m 2 z ) = 0.1207. Similarly as in the vector channel 
case we note a rather large difference between NLO and NNLO in the case of the MS 
scheme, which is significantly reduced if the preferred scheme is PMS. In NLO in 
the EC scheme we obtain 5?° = -0.0167 and A^L = 425 MeV (a s (m 2 T ) = 0.339). 



5 Discussion and conclusions 

Our discussion of the RS dependence of the QCD predictions and the fits may be 
summarized as follows: 

1. Changing the scheme from the MS scheme to the PMS scheme we obtain a 
reduction in the extracted value of a s {m^) by approximately 0.01 (a s (m|) is 
reduced by 0.001). Also, the difference between the NLO and NNLO results 
is much smaller in the PMS scheme than in the MS scheme. 
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2. Varying the scheme parameters r± and C2 in the region satisfying the condi- 
tion ( pOD with I = 2 we obtain an uncertainty in the extracted value of a s (m 2 ) 
of approximately 0.02 (uncertainty in a s (m 2 z ) is 0.002). As was argued in the 
text, this set of scheme parameters seems to be a minimal set that one should 
take into account in the discussion of the RS dependence. 

Our general conclusion is that the perturbative predictions for the QCD effects in 
the inclusive decay rates for the semileptonic r decays appear to be relatively precise, 
despite the rather low energy scale. It should be emphasized however, that in the 
discussion of the final precision of a s extracted from the r decays one should also 
take into account the approximate character of the SVZ expansion itself [^T], 53, |E5| . 



Discussions with R. Alemany, M. Beneke, M. Davier, A. Hocker, C. Maxwell and 
M. Neubert are gratefully acknowledged. 
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